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Origin of Unexpected Low Energy Structure in Photoelectron Spectra 
Induced by Mid-Infrared Strong Laser Fields 

Chengpu LiiQ an d Karen Z. HatsagortsyaqJ 
Max-Planck-Institut fur Kernphysik, Saupfercheckweg 1, D-69 117 Heidelberg, Germany 

' (Dated: July 30, 2010) 

Using a semiclassical model which incorporates tunneling and Coulomb field effects, the origin of the unex- 
pected low-energy structure (LES) in above-threshold ionization spectrum observed in recent experiments [C. I. 
Blaga et al, Nature Phys. 5, 335 (2009) and W. Quan et al, Phys. Rev. Lett. 103, 093001 (2009)] is identified. 
We show that the LES arises due to an interplay between multiple forward scattering of an ionized electron 
and the electron momentum disturbance by the Coulomb field immediately after the ionization. The multiple 
forward scattering is mainly responsible for the appearance of LES, while the initial disturbance mainly deter- 
mines the position of the LES peaks. The scaling laws for the LES parameters, such as the contrast ratio and the 
maximal energy, versus the laser intensity and wavelength are deduced. 

PACS numbers: 32.80.Rm; 32.80.Fb 



The physics of strong-field photoionization has been exten- 
sively studied in the last two decades and is assumed to be 
well-understood yj with the three-step model J2[] playing a 
fundamental role. The attosecond science has been developed 
based on the knowledge how to control the motion of elec- 
trons on the atomic scale which paved a way for new meth- 
ods for measuring molecular structure |3[|. The strong field 
approximation (SFA) |4] serves as an analytical tool to un- 
derstand the strong-field ionization phenomena which can de- 
scribe the rescattering via a perturbative approach [l|]. Re- 
cently, however, a number of strong-field ionization phenom- 
ena has been revealed which are caused by non-perturbative 
influence of the Coulomb field of the atomic core ||5|] : double- 
hump and interference structures Ja uD in the momentum 
distribution of photoelectrons near ionization threshold, frus- 
trated tunneling ionization |8|], and multiphoton assisted re- 
combination |9|]. With the advent and further improvement 
of intense femtosecond mid-infrared laser sources 11011 . the 
classical regime of strong-field ionization, when the Keldysh 
parameter y <C 1 and one expects the SFA to provide an ade- 
quate description, has been subjected to more attentive investi- 
gation. Here, y = y/I p /2U p , I p is the ionization potential, and 
Up the ponderomotive energy. The recent two experiments by 
C. Blaga et al. Ill 111 and W. Quan et al. 1121 on the photoion- 
ization of atoms and molecules in strong mid-infrared laser 
fields reveal a previously unexpected characteristic spike-like 
low-energy structure (LES) in the energy distribution of elec- 
trons emitted along the laser polarization direction, see Fig. 
1 (a). These observations manifest a striking contrast to the 
prediction of the SFA and point to a lack of complete under- 
standing of strong field physics. The numerical solutions of 
the Schrodinger equation within the single active electron ap- 
proximation can successfully reproduces the measured LES 
in the case of any atomic potential. However, the calculations 
using SFA or SFA with Coulomb corrections 11311 fail to de- 
scribe the LES. Varying the laser polarization from linear to 
circular, LES is significantly reduced. The latter indicates that 
forward rescattering [14] is playing an essential role in this 
process which has been pointed out by Blaga et al. 11 ill and 



Faisal 1 15]. However, the mechanism which creates the LES 
remains obscure. Many questions remain unanswered: How 
exactly does the LES arise? Why does it have a peaked struc- 
ture? Why is the effect of rescattering more pronounced in 
mid-infrared laser fields? 

In this Letter, we investigate in detail and identify the 
mechanism of LES. We employ the classical-trajectory Monte 
Carlo (CTMC) method with tunneling and the Coulomb field 
of the atomic core fully taken into account. In addition, we 
provide a qualitative theoretical estimation for the Coulomb 
field effects: initial Coulomb focusing (CF), multiple forward 
scattering and asymptotic CF. We quantify their relative role 
in the electron dynamics and conclude that 1) the behavior 
of the transverse (with respect to the laser polarization di- 
rection) momentum change of the electron due to Coulomb 
field effects with respect to the ionization phase is the key 
for understanding of the LES; 2) at mid-infrared wavelengths, 
multiple scattering of the ionized electron plays a decisive 
non-perturbative role. In particular, the transverse momen- 
tum change of the electron due to multiple scattering distorts 
the electron phase space to create peaks at low electron ener- 
gies, while the longitudinal momentum change due to initial 
CF shifts the peak energy to higher energies. We investigate 
and explain qualitatively the scaling of the LES parameters. 

In the 7< 1 regime, the electron oscillation amplitude 
in the laser field a = Eo/(0 2 exceeds the distance of the 
electron from the atomic core at the tunnel exit zo = Ip/Eo'. 
oc/zo ~ 2/y 2 ^> 1 and the transversal distance traveled by 
the electron during one laser period xq ~ 2^v^o/(2/ p ) 1 ' 4 0): 
Oc/xq ~ (l/27c)y/2I p /yco 7S> 1, where Eq and (0 are the laser 
field amplitude and frequency, respectively (atomic units are 
used throughout). Therefore, in this regime, the electron trav- 
els far from the core during its oscillation in the laser field, 
and the Coulomb field distorts the electron trajectory only at 
positions very close to the core. This happens immediately 
after ionization, corresponding to the initial CF, and when the 
electron revisits the atomic core and re-scatters. The num- 
ber of scatterings is large for mid-infrared wavelengths and 
low-energy photoelectrons 11611 : N s ~ a/xo ~ 10 at the pa- 
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FIG. 1: (color online). Photoelectron spectra (PES): (a) The ex- 
perimental result (squares) for a Xenon atom in a laser field with 
peak intensity Iq = 3.2 x 10 13 W/c/n 2 and wavelength X = 2.3pm 
UHl as well as the corresponding CTMC-simulation (circles), (b) 
CTMC-simulations for a hydrogen atom with A = 2pm, 1q = 9.0 x 
10 W/cm . exact (circles); with totally neglecting the Coulomb 
potential (NCP) (stars); NCP only for the electron longitudinal mo- 
menta P« (squares) and NCP only for the transverse momenta Pj_ 
(triangles). The high-energy limit of LES defined by the break in 
slope is indicated with an arrow. 



rameters in IT 1 1 M - The third Coulomb effect is the asymptotic 
CF when the electron momentum is disturbed by the Coulomb 
field after the laser pulse is switched off. This usually plays an 
important role for low-energy photoelectrons 1181. 1 171 1 1811 but is 
not essential for the LES as shown below. 

In our 3D CTMC simulation, the ionized electron wave 
packet is formed according to the ADK ionization rate 11911 
and further propagates classically. The electrons are born 
at the tunnel exit with the following conditions: (i) Along 
the laser polarization direction, the initial position Zi is de- 
rived from the effective potential theory 12011 and the initial 
momentum p t «, = 12111 : (ii) The transversal coordinates are 
Xj = yi = 0. The transverse momentum /?,-j_ follows the corre- 
sponding ADK distribution I19I1 . The transverse momentum 
components are p\ x = pi± sin (j) and pi y = pi± cos 0, where 
is the azimuthal angle randomly distributed within an interval 
of (Q,2n). The positions and momenta of electrons after the 
laser pulse are used to calculate the final asymptotic momenta 
111811 at the detector. Only electrons emitted along the laser 
polarization direction within an angle of ±2.5° are collected. 
The laser pulse is half-trapezoidal, constant for the first ten 
cycles and ramped off within the last three cycles. The elec- 
trons are launched within the first half cycle (of,- £ [0, n]). Our 
model provides an adequate description of the photoelectron 
spectrum (PES) and is qualitatively consistent with the exper- 
imental results as an example in Fig. \V[a) shows. 

Now we turn to the clarification of the physical mecha- 
nism behind the LES. For simplicity, we consider hydrogen 
atoms because the atomic structure is not essential for LES. 
Firstly, we investigate the changes in PES when some factors 
are neglected, see Fig. |T|(b). If the Coulomb potential is ne- 
glected after ionization, the LES disappears, as expected. If 
the Coulomb field effect is neglected only on the longitudi- 
nal momentum, the LES shifts towards lower energy and en- 
hances. However, if the Coulomb field effect is neglected only 
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FIG. 2: (color online), (a) The distribution of electrons in the LES in 
phase space (the initial transverse momentum versus the ionization 
phase <Pi = cot,) within the electron energy interval (0, 20) eV, with 
color coded probability. Irregular points are shown by triangles. The 
circles indicate the maximum probability for each phase. The laser 
field is maximal at <p,- = n/2. (b) PES with and without irregular 
points, (c) PES, the electrons born in the specified phase range are 
removed. The laser and atom parameters are the same as in Fig. 1(b). 



on the transverse momentum, the LES disappears completely. 
We can deduce that the change of the transverse momentum 
due to the Coulomb potential is the main source of LES at 
which we look more closely next. 

Secondly, we investigate the phase space distribution of 
electrons which contribute to the LES, the initial transverse 
momentum versus the birth phase shown in Fig. 12a): (i) 
The electron ionized near the peak of the laser field with a 
low transverse momentum has a low drift momentum, scatters 
many times by the atomic core with a low impact parameter 
and, as a result, shows chaotic behavior (the triangles in Fig. 
I2a)). However, the electrons with chaotic dynamics mostly 
contribute to the high-energy part of PES essential for the for- 
mation of a plateau in above-threshold ionization. The con- 
tribution of irregular points to the LES amounts only to a few 
percent 12211 . For this reason, when these points are discarded, 
the LES in PES shows only a little reduction in its lower en- 
ergy part, see Fig. |2b). Therefore, the chaotic dynamics is 
not the cause of the LES. (ii) A large number of electrons with 
a large initial transverse momentum (pi±) is concentrated af- 
ter the interaction in the phase space with a low transverse 
momentum. This is usually termed as CF 112.311 . The CF is 
largest near the laser peak and decreases with the increasing 
of the ionization phase, (iii) More importantly, the decrease of 
Pi± is not monotonic but shows a step-like slope change. We 
state that this slope change is responsible for the appearance of 
peaks in the LES. In fact, near the points of the slope change 
of the pi±((pi), the phase space of electrons contributing to the 
LES (per unit <p,-) has local maxima. This can be tested via ar- 
tificially discarding the electrons between certain characteris- 



tic phases, see Fig. EJc). The latter shows that the occurrence 
of the highest part of LES is, in fact, related to the electrons 
with ionization phases within an interval (1.64, 1.77) exactly 
corresponding to the slope change between B and D. The be- 
havior of p{± versus the ionization phase, in fact, reflects the 
one of the transverse momentum change due to the Coulomb 
field {8p±), see Fig. 2Ja) below. Then, one needs to explain 
8p± non-monotonic behavior with respect to 0;. 

Thirdly, we investigate the transverse and longitudinal mo- 
mentum disturbance of the electron due to the Coulomb field 
(5f>x and Spn). For one specific trajectory contributing to 
LES, 8p± and Sp\\ are shown in Fig. [2a). After the elec- 
tron's birth, Sp±\\ changes abruptly during a fractional part 
of the laser period which is known as initial CF [23]. Fur- 
ther momentum change takes place due to scattering when 
the electron passes the core at the minimum distance r (Fig. 
Ob)). We proceed to estimate the contributions of initial CF 
and multiple scattering, respectively, for electrons with the 
maximal probability for each phase (circles in Fig. 2(a)). 
(i) The transverse momentum change due to Coulomb po- 
tential V(r) at the sth scattering can be estimated as 8p± « 
f Vj_V(r(f ))dt ~ {p s /r^)8t s , where r s is the distance from the 
core at the scattering moment, p s is the one in the transversal 
plane, and 8t s is the scattering duration. When the electron 
velocity p s at scattering is large, 8t s ~ 2r s /p s . In the opposite 
case, 8t s ~ 2^2r s /\E((p s )\ is determined by the laser field 

E((p s ) at scattering. Accordingly, 8p± = —2p s /(r%p s ), if 

p] > r s \E(<p s )\, or 8p { l } = -2 3 / 2 p s /y/rj\E(<p s )\, otherwise. 

8p^ is sensitive to r s and p s . Their values for the different 
scattering events are shown in Figs. [3jc) and (d). As for the 




FIG. 3: (color online), (a) The transverse 5P± (solid) and longitudi- 
nal SP« (dashed) momentum changes due to the Coulomb field, (b) 
The distance r (solid) from the core and electric field (dashed) vs. the 
laser phase for one specific trajectory (<p; = 1 .72). (c) The position r s 
and (d) the momentum p s at the .sth scattering versus the ionization 
phase, (e) The estimated SPj_ due to different number of scatterings 
and initial CF. The laser and atom parameters are the same as in Fig. 
1(b). 
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FIG. 4: (color online). The estimation of (a) \SPj_\ and (b) SPu de- 
scribed in the text. The contributions of asymptotic CF the combined 
contribution of asymptotic and initial CF and the one of asymptotic 
CF and scattering are shown. The laser and atom parameters are the 
same as in Fig. 1(b). 
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formula of the Sp\ , p s in the expression for 8p ± should 
be substituted by the absolute value of the scattering coordi- 
nate |zj| along the laser polarization direction, (ii) The trans- 
verse momentum change due to initial CF can be estimated as 
8pf = -2p, ± / ± (£ , o,Co)|£'((p,)|/(2/p) 2 and the longitudinal 
one as Spl'*' = nf\\(E ,(o)\E((p i )\/(2I p ) ? '^, where the func- 
tions /i u(E(),(o) describe the deviation from the simple esti- 
mate of 118il . (iii) We estimate the asymptotic CF contribution 
comparing the asymptotic electron momentum with the one 
after switching off the laser pulse. 

In Fig. 3 (e), the relative contribution of multiple scattering 
(up to 4th scattering) and initial CF to 8p± is shown calcu- 
lated using the above formulas. At <p > 1.77 (line D), mainly 
1st scattering and initial CF plays a role. With a little shift 
of phase towards 1.68 (line C), the 2nd scattering enhances 
rapidly due to the increase of the scattering energy and the 1st 
scattering weakens. The significant contribution from the 2nd 
scattering results in the slope change between lines C and D 
in Fig. 2 (b). Decreasing the phase further, the 3rd scattering 
becomes comparable with initial CF, the 1st and 2nd scatter- 
ing. The competition between lst-3rd scattering and initial 
CF results in the flat slope between lines B and C (see Fig. 
2 (a)). The increasing contribution of the lst-3rd scattering 
as well as the emerging contribution from the 4th scattering 
are responsible for the further slope increase after the point B. 
Summing up our estimations for 8p± and 8pn due to scatter- 
ing, initial and asymptotic CF, the exact momentum change 
is reproduced, see Fig. 4. Neglecting scattering, initial or 
asymptotic CF, respectively, the contribution of each effect is 
quantified. The multiple scattering is crucial for the trans- 
verse momentum change. It determines the 8p±((pj) behavior 
which, in turn, determines the ionization phases correspond- 
ing to the LES peaks, and in this way the shape of LES. The 
initial CF plays a less important role for the transverse mo- 
mentum change, but a very significant one for the longitudinal 
momentum change. The latter shifts the LES peak, determin- 
ing the position of LES. These two points are consistent with 
the conclusions from Fig. 1(b). The asymptotic CF is only 
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FIG. 5: (color online). The dependence of En (circles) and R 
(squares) on (a) laser intensity Iq at A = 2fim, (b) wavelength A at 
/ = 9.0 x 10 13 W/cm 2 and (c) wavelength A at y = 0.534. The solid 
line is a fit to all of the LES data. The triangles and dashed lines are 
our estimations. In (d), (e) and (f) the electron distribution in phase 
space for different laser intensities and wavelengths are displayed. 



important at near zero energies and has little impact on LES. 

Finally, we investigate the dependence of LES on the laser 
intensity and wavelength, see Fig. 5. Two parameters charac- 
terize the LES, the first is Eh, the high-energy limit of LES 
1 1111 . as labeled in Fig. 1(b). The second is the contrast ra- 
tio (R), the ratio between the integral yield of photoelectrons 
in the energy interval (0,£7/) and the total ionization yield 
for electrons emitted along the laser polarization direction. R 
characterizes the relative height and Eh the width of the LES, 
respectively. Eh increases approximately linearly with inten- 
sity (see Figs. 5 (a)). Increasing the wavelength, Eh slightly 
deviates from the A 2 -law (i.e. the dependence on only 7) 
which is mainly due to the increasing initial CF contribution 
to the 5/7M . The R dependence on the laser parameters is more 
interesting. The R decreases monotonically with increasing 
intensity. With increasing wavelength, R first keeps constant 
and then decreases slowly. The Eh and R behavior can be ex- 
plained qualitatively by inspecting the phase space, see Figs. 
5 (d), (e). In Figs. 5(d) and (e), the curves in phase space ro- 
tate counter-clockwise resulting in the reduction of electrons 
within the LES. The reason is the competition between the 
initial CF and scattering. When fixing the Keldysh parame- 
ter (see Fig. 5(c), (f)), the curves in phase space shift paral- 
lel to larger phase which yields to the R increase with rising 
wavelength. At smaller wavelength, R, i.e. the LES visibility, 
becomes smaller. This partly explains why the LES has not 
been experimentally observed with near infrared lasers. An- 
other reason is the relative suppression of multiple scattering 
at near infrared wavelength due to quantum effects J22 \ . 

In conclusion, the experimental results have been success- 
fully reproduced via the semiclassical calculations and the ori- 
gin of LES is clarified. At mid-infrared wavelengths the multi- 
ple rescattering of ionized electron plays a significant nonper- 



turbative role which is the main factor creating the LES. The 
peaks in LES arise due to multiple scattering contributions to 
the transverse momentum. 

We gratefully acknowledge C. I. Blaga and L. F DiMauro 
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discussions. 
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